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Abstract 

Blind deconvolution is an ubiquitous non-linear inverse problem in applications like wireless communications 
and image processing. This problem is generally ill-posed since signal identifiability is a key concern, and there 
have been efforts to use sparse models for regularizing blind deconvolution to promote signal identifiability. Part I 
of this two-part paper establishes a measure theoretically tight characterization of the ambiguity space for blind 
deconvolution and unidentifiability of this inverse problem under unconstrained inputs. Part II of this paper analyzes 
the identifiability of the canonical-sparse blind deconvolution problem and establishes surprisingly strong negative 
results on the sparsity-ambiguity trade-off scaling laws. Specifically, the ill-posedness of canonical-sparse blind 
deconvolution is quantified by exemplifying the dimension of the unidentifiable signal sets. An important conclusion 
of the paper is that canonical sparsity occurring naturally in applications is insufficient and that coding is necessary 
for signal identifiability in blind deconvolution. The methods developed herein are applied to a second-hop channel 
estimation problem to show that a family of subspace coded signals (including repetition coding and geometrically 
decaying signals) are unidentifiable under blind deconvolution. 

Index Terms 

Identifiability, rank one matrix recovery, blind deconvolution, parametric representation, rank two null space 


1. Introduction 

B lind deconvolution is a challenging open problem in signal processing that manifests itself in a host of 
applications like image de-blurring, blind system identification, and blind equalization. The difficulty of blind 
deconvolution stems from its ill-posedness in the sense of Hadamard [3], [4] as well as non-linearity of the problem. 
In practice, additional application specific prior information on the unknown signals is necessary to render this 
inverse problem better behaved. Sparsity based models have been used extensively in the past decade to capture 
hidden signal structures in many applications of interest. While there have been a few attempts at exploiting sparsity 
based priors for blind deconvolution type of problems [5]-[10], strong theoretical guarantees on signal identifiability 
for such sparse models are largely missing with the exception of [5] (which proves an interesting positive result under 
idealized assumptions). Exploring more realistic sparsity assumptions pertaining to multi-hop channel estimation, 
we prove some surprising negative results for signal identifiability in the present article. Part I of this paper [II] 
was devoted to the characterization of the ambiguity space associated with blind deconvolution and the proof of a 
strong unidentifiability result. Part II of the paper utilizes the results in Part I to characterize the ambiguity space for 
sparsity constrained blind deconvolution problems. In particular, canonical-sparse blind deconvolution is considered 
and a trade-off is derived between the level of sparsity in the unknown signals and the dimension of the ambiguity 
space. Unfortunately, the derived scaling laws turn out to be unfavorable when compared to analogous results 
in compressed sensing [12]. Furthermore, an application of the methodology to a theoretical abstraction of the 
second-hop channel estimation problem, yields an unfavorable sparsity-ambiguity trade-off for a family of subspace 
coded signals that include repetition coding and geometrically decaying signals. 
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A. Related Work 

We briefly reviewed prior research on the multi-channel version of blind deconvolution [13]-[15] in Part I (see [11] 
for details). We briefly recap prior art discussed in Part I of the paper that pertains specifically to blind deconvolution 
in single-input-single-output (SISO) systems, before mentioning more related research. 

For SISO blind deconvolution {i.e. blind channel estimation), [16] showed that input signal realizations with 
no inter-symbol-interference (ISI) could be uniquely recovered up to scalar multiplicative ambiguities. In contrast, 
we showed in Part I of the paper that if ISI were allowed, blind deconvolution is ill-posed to the point of almost 
every input signal pair being unidentifiable. In the papers [5], [17], blind deconvolution was cast as a rank one 
matrix recovery problem for the development of convex programming based computational heuristics and a proof 
of correctness was provided assuming a random subspace precoding with sufficient redundancy. We adopted this 
framework in our earlier works [18]-[20] on an information theoretic characterization of identifiability of the 
solution to general bilinear inverse problems and showed good predictive performance of the theory for blind 
deconvolution with Gaussian/Rademacher inputs. Whereas [20] derived conditions for identifiability (analogous 
to achievability results in information theory) and developed an abstract theory of regularized bilinear inverse 
problems, making assumptions no more specific than non-convex cone constrained priors, the present paper makes a 
significant departure in the sense that more explicit sparsity and subspace priors are considered that are motivated by 
applications in multi-hop channel estimation [21]-[23] and the focus is on characterizing the inherent unidentifiability 
of blind deconvolution (akin to converse/impossibility results in information theory). We point out that beside 
multi-hop communications, convolution with sparse vectors is an intuitive modeling assumption in other applications 
as well, e.g. high resolution astronomical imaging [9], modeling echoing audio signals [24], ultrasound imaging [25], 
and sporadic communications in 5G systems [26]. 

A promising identifiability analysis was proposed in [6], leveraging results from [27] on matrix factorization for 
sparse dictionary learning using the ^i-norm and Iq quasi-norm for 0<Q< 1. Their approach and formulation differ 
from ours in two important aspects. Firstly, we are interested in single-in-single-out (SISO) systems whereas [6] deals 
with single-in-multiple-out (SIMO) systems. Secondly, [6] analyzes identifiability as a local optimum to a non-convex 
ii (or iq for 0 < q < 1) optimization and hence is heavily dependent on the algorithmic formulation, whereas we 
consider the solution as the local/global optimum to the Iq optimization problem and our impossibility results are 
information theoretic in nature, implying that they hold regardless of algorithmic formulation. We emphasize that 
the constrained ii optimization formulation in [6] is non-convex and therefore does not imply existence of provably 
correct and efficient recovery algorithms, despite identifiability of the channel. Although it would be interesting to 
try and extend their approach to SISO systems and compare with our results, this is non-trivial and beyond the 
scope of the present paper. 

An inverse problem closely related to blind deconvolution is the Fourier phase retrieval problem [28]-[30] 
where a signal has to be reconstructed from its autocorrelation function. This is clearly a special case of the blind 
deconvolution problem with much fewer degrees of freedom and allows identifiability and tractable recovery with a 
sparsity prior on the signal [28]. A second important difference is that after lifting [31], the Fourier phase retrieval 
problem has one linear constraint involving a positive semidefinite matrix. This characteristic is known to be helpful 
in the conditioning of the inverse problem and in the development of recovery algorithms [32]. While the blind 
deconvolution problem does not enjoy the same advantage, this approach seems to be a good avenue to explore if 
additional constraints are allowed. This is a future direction of research. 

As in Part I, we rely on the parameter counting heuristic to quantify ‘dimension’ of a set/space. The parameter 
counting heuristic can be made rigorous using the framework of Hausdorff dimension [33], but we shall not bridge 
this gap in the present paper owing to space limitations. 

B. Contributions and Organization 

In this work, we quantify unidentifiability for certain families of noiseless sparse blind deconvolution problems 
under a non-asymptotic and non-statistical setup. Specifically, given model orders m^n ^ Z+, we investigate the 
trade-off between the level of sparsity in the unknown signal pair G /C C x and the dimension 

of the unidentifiable subset of signals in 1C. While Part I of the paper focused on developing a description of the 
ambiguity space. Part II studies the interaction of this ambiguity space with canonical-sparse vectors and further 
extends the techniques to consider repetition coded and geometrically decaying vectors. To motivate our signal 


3 


choices, the specific application we consider is multi-hop sparse channel estimation for relay assisted communication 
described in Section IV. Like Part I, our focus is on an algorithm independent identifiability analysis and hence we 
shall not examine efficient/polynomial-time algorithms, but rather show information theoretic impossibility results. 
Our approach leads to the following novelties for Part II (refer to [ ] for contributions of Part I of the paper). 

1) We show that sparsity in the canonical basis is not sufficient to ensure identifiability, even in the presence of 
perfect model order information, and we construct non-zero dimensional unidentifiable subsets of the domain 
/C for any given support set of (or y^) as evidence to quantify the sparsity-ambiguity trade-off. This is the 
content of Theorem 2 in Section III-B. 

2) We consider a theoretical abstraction to a multi-hop channel estimation problem and extend our unidentifiability 
results to this setting. Specifically, we show that other types of side information like repetition coding or 
geometric decay for the unknown vectors are still insufficient for identifiability, and exhibit an ambiguity 
trade-off analogous to the case of canonical-sparsity. This is the content of Theorem 3 and Corollaries 3 and 6 
in Section IV. 

The rest of the paper is organized as follows. Section II-B presents constructive numerical examples illustrating 
the nature of unidentifiable inputs for blind deconvolution. Section II-A recaps the system model, the notion of 
identifiability, and the lifted reformulation of the blind deconvolution problem from Part I of the paper. Section III 
presents the key unidentifiability results for canonical-sparse blind deconvolution and contrasts them against results 
from Part I of the paper. Section IV extends the techniques and results in Section III to more general families of 
subspace based structural priors, including repetition coding and geometrically decaying canonical representations. 
Section V concludes the paper. Detailed proofs of the major results in the paper appear in Appendices A-C. 


C. Notational Conventions 

We shall closely follow the notation used in Part I of the paper which is recapped below. All vectors are assumed 
to be column vectors unless stated otherwise and denoted by lowercase boldface alphabets {e.g. a). Matrices are 
denoted by uppercase boldface alphabets {e.g. A). The MATLAB® indexing rules are used to denote parts of 
a vector {e.g. a(4 : 6) denotes the sub-vector of a formed by the 4^^, 5^^ and 6^^ elements of a). The all zero 
vector/matrix (respectively all one vector) shall be denoted by 0 (respectively 1) and its dimension would be clear 
from the usage context. For vectors and/or matrices, (•)^ and rank( ) respectively return the transpose and rank 
of their argument, whenever applicable. Special sets are denoted by uppercase blackboard bold font {e.g. M for 
real numbers) and a subscripted ‘+’ sign would denote the non-negative subset whenever applicable {e.g. R+ for 
non-negative real numbers). Other sets are denoted by uppercase calligraphic font {e.g. S). For any set S, |5| shall 
denote its cardinality. Linear operators on matrices are denoted by uppercase script font {e.g. 

To avoid unnecessarily heavy notation, we shall adopt the following convention: The scope of both vector variables 
(like X, y, u, v, etc.) as well as matrix variables (like X, Y, etc.) are restricted to individual theorems and/or 
proofs. Their meanings are allowed to differ across theorems and proofs (and even across disjoint subparts of 
the same proof when there is no risk of confusion), thus facilitating the reuse of variable names across different 
theorems and avoiding heavy notation. 


IT Background 

In this section, we review the system model from Part I of the paper and overview a new result from Part II of 
the paper with a numerical example. The results in Part II of the paper continue to use the same system model, but 
with a different feasible set of signals. 

A. System Model 

We consider the noiseless linear convolution observation model 


z — xi^y^ 


( 1 ) 
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where ik: x R’^ ^ ]^rn+n i denotes the linear convolution map, {x,y) e )C C R^ x R’^ denotes the pair of 

input signals from a restricted domain /C, and z G is the vector of observations given by 


^ x{j)y{l + l-j), l<l<n, 


z{l) = < / ,} . 


( 2 ) 


x{j)y{l + 1 — j), 1 < / — n < m — 1. 


The corresponding inverse problem of constrained blind linear deconvolution is to find the vector pair (x,y) G /C 
given the noiseless observation z and is symbolically represented by the feasibility problem 


find {x,y) 

subject to X X y = z, (Pi) 

{x,y) G JC. 


Assuming that the model orders m and n, respectively, of vectors x and y are fixed and known a priori, we are 
concerned with whether Problem (Pi) admits a unique solution. Since bilinearity of the convolution operator implies 
the inherent scaling ambiguity in Problem (Pi) given by 

xi^y — axi^—y^ Va 0, (3) 

cr 

the question of uniqueness only makes sense modulo such scalar multiplicative factors. This leads us to the definition 
of identifiability in Part I that we restate below. 


Definition 1 (Identifiability). A vector pair (x,y) G /C C R’^ x R’^ is identifiable within 1C with respect to the 
linear convolution map x, if\l{x'^ y') G 1C satisfying xi^y — x' i^y', 3a 0 such that (x', y') = (ax, ^y\ 

It is easy to see that Definition 1 induces an equivalence structure on the set of identifiable pairs in 1C. For 
future reference, we define the equivalence relation IdRi /C x /C ^ {0,1} as follows. Given any (x, y), (x', y') G 1C, 
IdR((x, y), (x', y')) = 1 if and only if 3a 0 such that (x',^) = (ax, ^y). It is straightforward to check 
that IdR(-, •) is indeed an equivalence relation. Let /C/Mr denote the set of equivalence classes of K induced 
by IdR(-, •), and for any (x,y) G 1C let [(x,y)] G /C/IdR denote the equivalence class containing (x,y). Then 
Definition 1 amounts to declaring a vector pair (x,y) G /C as identifiable if and only if every (x',y') G 1C with 
[(x', y')] [(x, y)] satisfies xi^y ^x' i^y'. 

Using the lifting technique from optimization [31], Problem (Pi) can be reformulated as a rank minimization 
problem subject to linear equality constraints [17], [18] 

minimize rank ( W) 

w 

subject to y{W) = z, (P 2 ) 

TU G W, 


where W C is any set satisfying 

>Vp|{W e I rank(W) < 1} = {xy'^ \ {x,y) E JC}, (4) 

and ^—)■ is the unique linear operator (henceforth referred to as the lifted linear convolution 

operator) satisfying 

^{xy'^) ^xi<y, V(®,y) € M"* X M** (5) 

and admitting an explicit closed form specification (see Part I [11], Remark 1). Problems (Pi) and (P 2 ) are equivalent 
for z 0 (see [20] for details). Hence, owing to the analytical simplicity of Problem (P 2 ), identifiability of the 
solution to Problem (Pi) is derived using the properties of Problem (P 2 ). By construction, the optimal solution to 
Problem (P 2 ) is a rank one matrix Wopt and its singular value decomposition Wopt = cTopt'^^opt'^^t yields a solution 
(^7 ?/)opt ~ (v^^opt'^opt? sj ^opt'^opt) fo Problem (Pi). 
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B. An Overview of the Results 

Let us interpret the role of /C C x for the discrete-time blind linear deconvolution problem (Pi). Suppose 
that G /C is the ground truth resulting in the observation z = The set 1C captures application 

specific constraints on the signal pair (in our case, sparsity constraints on the true solution (x*,y*)) by 

restricting the feasible set of signal choices in Problem (Pi). Without such a restriction, the solution to Problem (Pi) 
would most likely be unidentifiable (infact. Theorem 2 in Part I shows that almost every ground truth {x^^yf) 
w.rt. the Lebesgue measure is unidentifiable in the absence of constraints). In the sequel, we solely consider fC 
to be a separable cone in x M^, Le. there exist sets Vi C and V2 C MA such that )C = Vi x V2 and for 
every cr > 0, {x^y) G /C implies that (ax^ay) G 1C. Such a separability of 1C is motivated by the observation that 
in many applications of interest x^ and y^ are unrelated (like x^ and y^ may respectively represent the unknown 
source and the unknown channel in a blind channel estimation problem). 

Theorem 3 in Section IV-B is the most general result of this part of the paper. However, we shall not state it here 
since it requires some technical definitions that would appear unmotivated at this stage. Instead, we shall succinctly 
overview the results pertaining to the special case of canonical-sparse blind deconvolution. We revisit the numerical 
example from Part I of the paper and examine its identifiability (according to Definition 1) under discrete-time blind 
linear deconvolution with canonical-sparsity constraints. In preparation for the result, we first define a parameterized 
family of canonical-sparse cones for an arbitrary integer > 3 and any index set A C {2, 3,..., (i — 1} as 


K:o{A,d) ^ 


iu(l) ^ Q,w{d) ^ 0,iu(A) = o|, 


( 6 ) 


i.e. A denotes the set of indices that are zero across all vectors in /Co (A, d). Secondly, let us denote the Minkowski 
sum of the sets { — 1} and A, by the shorthand notation A — 1, defined as 


A - 1 = {-1} + A ^ {j - 1 I i € A}. 


(7) 


Jumping ahead to provide a concrete example, the following unidentifiability result can be obtained as a corollary 
to Theorem 2 in Section III-B. 


Corollary 1. Let m > 5 and n > 2 be arbitrary integers and 0 ^ A C {3,4,..., m — 2} denote a set of indices. 
Let 1C = /Co(A,m) x be the structured feasible set in Problem (Pi) and define p = |A|J(A — 1)|. Then there 
exists a set C /C/ Mr of dimension {m + n — p — 1) such that every signal pair (x, y) G is unidentifiable by 
Definition 1. 


Corollary 1 cannot be checked in its full generality by considering examples, since it is a statement about an 
uncountably infinite number of vectors. However, we will re-analyze the numerical example in Part I of the paper 
(within canonical-sparse feasible sets) to capture the ideas behind Corollary 1. As in Part I, consider the vectors 

= (1,0, 1 , 0 , 0 , 0 , 0 , 0 ,1,0,1)^, = (1,0,0,0,1,0,0)^, (8a) 

(8b) 


resulting in 


a :2 = ( 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , 0 , 0 )"^, ?/2 = ( 1 , 0 , 1 , 0 , 1 , 0 , 1 )"^, 
xi-kyi^ X2-ky2^ ( 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 1 , 0 , 0 )^ 


(9) 


with xi and X 2 being non-collinear. Clearly, the pairs (xi,yi) and (x 2 ,y 2 ) are unidentifiable within the domain 
X M^. Setting = 11 and A = {4, 5,6, 7, 8} in (6) gives xi, X 2 G /Co(A, 11) and therefore (xi, yi), (x 2 , 2 / 2 ) ^ 

/Co (A, 11) xMJ are still unidentifiable within /C = /Co (A, 11) x . 

Next, we show that the rotational ambiguity over x (discussed in Part I) stays valid over /Co(A, 11) x R^ 
and leads to an uncountable set of unidentifiable pairs in 1C. Let zq = xix yi = X 2 X y 2 denote the common 
convolved output in (9) and consider the parameterized vectors 

x'l = xi cos 9 — X2 sin y[ = yi sin f — y2 cos 0, (10a) 

X 2 = xi cos — X 2 sin 0, ^2 — yi sin 9 — y 2 cos 0, (10b) 


where 0 ^ 0 are the parameters, and {xi, X 2 , yi, ^ 2 } acts as the set of seed vectors for the above transformation. 
Clearly, 9 ^ f and non-collinearity of xi and X 2 imply that x'^ and X 2 are linearly independent. A simple algebraic 
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manipulation reveals that 

^ 2/i = ^2 ^ 2/2 — ^0 sin(0 + 0) — a^2 ^ 2/1 sin 0 sin 0 — xi ★ ^2 cos 9 cos 0, (11) 

rendering both and ( 2 ^ 252 / 2 ) unidentifiable within x Since xi,X 2 ^ ^o(A, 11), (10) implies that 

X 2 G /Co(A, 11), thus rendering ^ 2 ) ^ A;o(A, 11) x unidentifiable within 1C = /Co(A, 11) x 

Since (6^,0) G [0,7r)^ describes a two dimensional parameter space, (10) and (11) imply that the unidentifiable 
subset of /Co(A, 11) xM7 is at least two dimensional and therefore hints towards Corollary 1. Thus, having a sparse 
support does not help with identifiability here. 

III. Unidentifiability under Canonical Sparsity 

We shall use identifiability in the sense of Definition 1. In Section III-A, we recap the partially parametric 
characterization of the ambiguity space of unconstrained blind deconvolution (Lemma 1), the pathological cases for 
unidentifiability, and a rotation based representation cum decomposition result (Lemma 2) discussed in Part I of 
the paper. In Section III-B, we state our main unidentifiability result for canonical-sparse blind deconvolution as 
Theorem 2 and contrast it with the almost everywhere unidentifiability result for non-sparse blind deconvolution 
from Part I of the paper. Section III-C states a stronger result for the feasible set in Corollary 1. Throughout this 
section, we assume that /C represents the (not necessarily convex) feasible cone in Problem (Pi), i.e. V(x,y) G 1C 
one has (dx, ay) G 1C for every a ^ 0, but 1C is allowed to change from theorem to theorem. 


A. Representation Lemmas and Pathological Cases 

Let y{') denote the lifted linear convolution operator as described in Section II-A. We denote the rank-fc null 
space of by k) and define it as 

k) = {Qe I rank(Q) < fc, ^(Q) = O}. (12) 

We note that the rank one null space of y{-) is trivial, i.e. 1) = {0}, and this property of linear convolution 

is implicitly used to prove equivalence of Problems (Pi) and (P 2 ) for z ^ 0 in [20]. A/’(^, 1) = {0} follows from 
interpreting convolution as polynomial multiplication, since the product of two real polynomials is identically zero 
if and only if at least one of them is identically zero. The following lemma (borrowed from Part I of this paper) 
describes a subset of 2), the rank two null space of y{-), and is used in the proofs of the results in the 

sequel. 


Lemma 1 (from [11]). Let m,n > 2 and Q G admit a factorization of the form 

0 

0 J ’ 

for some v G and u G Then Q G 2). 




u 0 
0 —u 


(13) 


We refer the reader to Remark 2 in [11] for correct parsing of the symbolic shorthand on the r.h.s. of (13). We 
shall also need the following non-linear re-parameterization cum decomposition result (borrowed from Part I of this 
paper) to serve as a building block for constructing adversarial instances of input signals for which deconvolution 
fails the identifiability test. We point out the symbolic connection to the transformation in (10) and the representation 
in (13). 


Lemma 2 (from [11]). Let d > 2 be an arbitrary integer and w G [w' G R^ | w'{l) ^ O^wfd) ^ O} be an 
arbitrary vector. The quotient set Qr^{wy) defined as 




|(m*,7) G R^ ^ X [0,27r) 



0 

-w^ 


cos 7 
sin 7 


(14) 


is finite (possibly empty) with cardinality at most {2d — 2). Ifd is an even integer then Q^(m, d) is non-empty. 


We briefly recap the pathological cases discussed in Part I, Section IV-A that we shall exclude from consideration 
in rest of the paper. This serves to remove the most straightforward delay ambiguities at an intuitive level, as 
well as to simplify the arguments in our proofs by excluding special cases that are best considered separately. 
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0 2 4 6 8 10 12 14 

index, I 


Fig. 1. An arbitrary vector s G /Co(A, d) with A = {3,4, 7, 8,9,12} and d = 14. Every s G /Co(A, d) is zero on the index set A (indicated 
by blue dots). Heights of the black dashed stems, indicating values on A^, can vary across different vectors in /Co(A,d). 


It was shown in Part I that for identifiability of (x*, y^) G x within x it is necessary that both 
^ {0} and {x*(l), ^ {0} must be true. In the sequel, we shall consider the stronger 
restriction 0 0 {x*(l), x*(m), y*(l), to automatically eliminate the pathological cases. Implicitly, this is 
also the reason for requiring w{l) ^ 0 and w{d) ^ 0 for any vector ru G in the premise of Lemma 2 as well as 
in the definition of the canonical-sparse domain /Co(A,(i) in (6). 


B. Canonical-sparse Blind Deconvolution 

To be consistent with the notation for canonical-sparse cones in (6), we denote the set of unconstrained non- 
pathological d dimensional vectors by 


/C(0,d) ^ 


^ 0,w{d) ^ o|. 


(15) 


For subsequent comparison we recall below, the unidentifiability result from Part I of the paper adapted for the 
domain /C(0, (i). 


Theorem 1 (adapted from [11]). Let m^n > 4 be even integers and 1C — /C(0,m) x W^. For any x G /C(0,m), 
(x,y) G fC is unidentifiable almost everywhere w.rt. any measure over y that is absolutely continuous w.r.t. the n 
dimensional Lebesgue measure. 


In the presence of a sparsity prior on y G Theorem 1 does not apply anymore, since a sparsity prior is 
necessarily generated from a measure that is not absolutely continuous w.r.t. the n dimensional Lebesgue measure. 
Assuming that the sparsity prior is w.r.t. the canonical basis, we prove the following unidentifiability result (recall 
that /C/Mr denotes the set of all equivalence classes induced by the equivalence relation IdR(-, •) on /C). 

Theorem 2. Let m^n > 5 be arbitrary integers. For any given index sets 0 ^ Ai C {3,4,...,m — 2} and 
0 7 ^ A 2 C {3,4,..., n — 2}, let 1C — /Co(Ai, m) x /Co(A 2 , n) and define pj = \Aj U(Aj — 1)| for j G {1, 2}. Then 
there exists a set C /C/ Mr of dimension {m-\- n — 1 — pi — P 2 ) such that every (x, y) G is unidentifiable. 

Proof: Appendix A. ■ 

As a visualization aid, an arbitrary vector s in the canonical-sparse cone /Co(A, d) = /Co({3,4, 7, 8 ,9, 12 }, 14) is 
shown in Fig. 1. A couple of comments about the premise of Theorem 2 are in order. Let (x*, y*) G /C denote an 
unidentifiable pair in the feasible set of Theorem 2. 

1 ) Our proof technique hinges on an adversarial construction of x*(Ai |J(Ai — 1)) and y*(A 2 U (^2 — !))• By 

feasibility of {x^{l)^x^{m)^y^{l)^y^{n)} ^ 0 needs to be satisfied, which might conflict with our 

adversarial construction if {1, m) n(Ai U(^i “ 1)) 7 ^ 0 0(^2 U (^2 “ 1)) 7 ^ 0 holds. Clearly, both 

of these scenarios are rendered impossible, if we insist on Ai C {3,4,..., m — 2} and A 2 C {3,4,..., n — 2}. 

2) We insist on Ai ^ 0 and A 2 7 ^ 0 to have a strictly non-trivial realization of the sparse blind deconvolution 
problem, i.e. an instance which violates some assumption of Theorem 1 other than model orders m and n 
being even. It is easy to see that if A = 0 in ( 6 ) then /Co(A, d) = /C(0, d) admits a non-zero d dimensional 
Lebesgue measure. Hence, if either Ai or A 2 is empty in Theorem 2, then the sparse blind deconvolution 













problem instance so generated will fall under the purview of Theorem 1 (assuming even m and n). Technically 
speaking, Theorem 2 is valid even if Ai or A 2 is empty (provided that m^n > 4), however its implications are 
weaker than that of Theorem 1 when m and n are even. Furthermore, note that the requirements Ai 0 and 
A 2 0 respectively imply m > 5 and n > 5. 

We note that the assumptions of Theorem 2 imply pi < {m — 3) and P 2 < — 3) so that the unidentifiable subset 

of 1C/ ITr is at least 5 dimensional. In particular, the unidentifiable set in Theorem 2 is always non-trivial. We also 
note that the canonical-sparse feasible domain of Theorem 2 is far more structured than the feasible set of Theorem 1 
and hence, the set of all unidentifiable inputs in Theorem 2 is much smaller. Nonetheless, the canonical-sparsity 
structure is not strong enough to guarantee identifiability for all canonical-sparse vectors. Additionally, note that 
Ai and A 2 denote sets of zero indices, so that larger cardinality of Ai or A 2 implies a sparser problem instance. 
Furthermore, we note that Corollary 1 straightforwardly follows from Theorem 2 by setting Ai = A and A 2 = 0. 
However, the result in Corollary 1 is not tight for n > 4 being an even integer. We remedy this shortcoming in 
Section III-C by directly specializing Theorem 1 to /C = /Co(A,m) x (the feasible set of Corollary 1) to yield 
Corollary 2. 

Exploiting the equivalence between bilinear inverse problems and rank one matrix recovery problems [20] (in 
the current context, equivalence between Problems (Pi) and (P 2 )) the result of Theorem 2 can be interpreted as 
evidence of the null space of the convolution operator admitting a large number of simultaneously canonical-sparse 
and low-rank matrices. For rank one matrix completion problems [34], it is relatively straightforward to see that 
a random sampling operator on a sparse rank one matrix will return zeros on most samples, thus rendering it 
impossible to distinguish the rank one matrix in question from the all zero matrix. However, the same observation 
is not at all straightforward for a rank one matrix recovery problem [35] when the sampling operator is fixed to 
the lifted linear convolution operator. Theorem 2 asserts that this is indeed true and the lifted linear convolution 
operator y{-) admits a large number of non-zero canonical-sparse matrices within its rank two null space. 

It is important to note that Theorems 1 and 2 are of different flavors and are not comparable since they make 
different assumptions on the feasible domain. In particular, neither theorem universally implies the other even if we 
consider special cases for each of them. To make this point more explicit, we make the following observations. 

1) Theorem 1 asserts almost everywhere unidentifiability within the feasible set /C(0,m) x M^, but this does not 
imply the conclusions of Theorem 2, since the feasible set in the latter theorem is /Co(Ai,m) x /Co(A 2 ,n), 
which is a measure zero set w.r.t. the (m + n) dimensional Lebesgue measure associated with the Cartesian 
product space /C(0,m) x W^. 

2) In Theorem 2, even if we assume the model orders m and n to be even and the zero index sets Ai and A 2 
to be empty (so that the feasible set /Co(Ai,m) x /Co(A 2 ,n) is equal to /C(0,m) x almost everywhere 
w.r.t. the (m + n) dimensional Lebesgue measure), we can only draw the conclusion that there exists a (m + n) 
dimensional unidentifiable subset (call it ^unio) of the feasible set. This is clearly insufficient to show almost 
everywhere unidentifiability (as claimed by Theorem 1) since the complement of ^uniD could be a (m + n) 
dimensional set as well, admitting non-zero (m + n) dimensional Lebesgue measure. 

3) If Ai = A 2 = 0 is considered with even model orders m,n > 4 so that the feasible sets in Theorems I 
and 2 are equal almost everywhere w.r.t. the (m + n) dimensional Lebesgue measure, then the conclusion of 
Theorem I is stronger, since almost everywhere unidentifiability over /C(0,m) x automatically implies the 
existence of a (m + n) dimensional unidentifiable subset. 

C. A Mixed Extension 

It is possible to fuse the ideas from Theorems 1 and 2 to strengthen the conclusion of Corollary 1 on the domain 
^ = ^o(A, m) X i.e. on a domain formed by the Cartesian product of a canonical-sparse and a non-sparse cone. 

Corollary 2. Let m > 5 be an arbitrary integer and n > 4 be an even integer. For any given index set 9 A C 
{3,4,..., m — 2}, let 1C — /Co(A, m) x and define p = \A IJ(A — 1)|. Then there exists a set C /Co(A, m) 
of dimension {m — p) such that Vx E Q', (x^y) E 1C is unidentifiable almost everywhere w.r.t. any measure over y 
that is absolutely continuous w.r.t. the n dimensional Lebesgue measure. 


Proof: Appendix B. 
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Fig. 2. Two hop relay assisted communication link topology between the source S and the destination D with k intermediate relays Ri, R 2 , 
Rk. The SISO channel is scaled by g{lj) and delayed by Ij units (represented as g(lj)D~^^h) for l<j<k. 


Note that Corollary 2 is a non-trivial extension of Theorem 1, since Definition 1 defines identifiability of an 
input signal pair within some feasible set 1C and this set is different for Theorem 1 and Corollary 2. In particular, if 
G /C is an unidentifiable pair in the feasible set of Corollary 2, then the proof of the result involves the 
construction of a candidate adversarial input G /Co (A, m) as in Theorem 2 as well as the vector G as in 

Theorem 1. 

Through the unidentifiability results in this section, we have attempted to quantify the ill-posedness of blind 
deconvolution under canonical-sparsity constraints, and illustrate the geometric reasons behind this unfavorable 
sparsity-ambiguity trade-off behavior. The rule of thumb for determining the sparsity-ambiguity trade-off for 
canonical-sparse blind deconvolution is well illustrated by comparing the conclusions of Theorem 1 and Corollary 2. 
Roughly speaking, if the set of zero indices for the signal vector x G has cardinality |A| then the dimension of 
the set of unidentifiable choices for x is reduced from m in Theorem 1 to somewhere between (m — |A| — 1) and 
(m — 2|A|) in Corollary 2. The main message is that if an application exhibits a bilinear observation model of linear 
convolution and the additional application specific structure is that of canonical sparsity of the unknown variables, 
then it is necessary to drastically revise the system design specifications for any hope of signal identifiability. Such 
revision may incorporate some form of randomized precoding of the unknowns, so that the effective bilinear operator 
governing the observation model looks substantially different from the convolution operator {e.g. the Gaussian 
random precoding used in [5]). Alternatively, sparsity in non-canonical bases could also be helpful {e.g. Rademacher 
random vector signal model used in [20] and the constant modulus based problem instances in [14]). 

IV. Unidentifiability Results with Coding 

Structural information about unknown input signals is almost exclusively a result of specific application dependant 
data or system deployment architectures. In this section, we consider an abstraction for a problem in multi-hop 
channel estimation [21], [22] and analyze other subspace based structural priors relevant to this application. Studying 
constrained adversarial realizations of these abstractions has utility far beyond communication systems. For example, 
when communication subsystems are parts of a larger complex system involving sensing and navigation for 
autonomous vehicles [36]-[40] then non-linear inter-system interactions may become unavoidable. In this section, 
we extend our unidentifiability results from Section III to simple forms of coding across the transmitters for the 
second-hop channel estimation problem. Section IV-A considers the somewhat idealized, but simpler to interpret, 
scenario of repetition coding across transmitters and states unidentifiability results for both unstructured as well as 
canonical-sparse channels. Rather than explicitly proving the results in Section IV-A, we note that they are important 
special cases of the more general unidentifiability result in Section IV-B for certain families of partially cooperative 
codes (Theorem 3) that we prove in Appendix C. Finally, Corollary 6 states yet another important special case of 
Theorem 3 pertaining to geometrically decaying subspace sparse signals, previously considered in [2]. In the spirit of 
using commonly employed notation, we shall denote the unknown pair {x^y) by the pair (^, h) with g pertaining 
to a topological configuration of the network and h representing the channel impulse response as described next. 

In a previous paper [21], we developed a structured channel model based on the multichannel approximation 
for the second hop in a relay assisted communication topology shown in Fig. 2. Specifically, we consider k 
relays Ri, R 2 , ..., R^ and let h denote the propagation delay and power adjusted SISO channel impulse 
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response (CIR) common to each relay-destination Rj D channel, j G {1, 2,..., k}. Thus, the effective CIR at 
the destination is formed by a superposition of k distinct delayed scalar multiples of the vector h. If the maximum 
propagation delay of Rk relative to Ri is upper bounded by m, then there exists a vector g G and an index subset 
{/i, / 2 , • • •, 4} ^ {I 5 2,..., m} (with li = 1) such that the propagation delay adjusted SISO CIR for the Rj ^ D 
channel is g{lj)h and that g{l) = 0, V/ 0 {hj 2 ^ • • •, h}- With these definitions and letting G {0, 
denote the Toeplitz matrix representation for delay by I units VI < / < m, the effective CIR at the destination is 
given by the linear convolution z = Yl^=i = g i^h ^ channel estimation problem at 

the destination is to recover the vector pair (^, h) from the observation of z. Since the vectors g and h are not 
arbitrary, but have physical interpretations leading to structural restrictions, we can incorporate this knowledge by a 
constraint of the form of (g^h) G /C for some set JC C x W^. For example, if h represents a SISO CIR for 
underwater acoustic communication [41], [42] (or, more generally, for wide-band communication [43]), then h has 
been observed to exhibit canonical-sparsity. Now, associating the pair (^, h) with a feasible solution to Problem (Pi), 
the success of the channel estimation problem at the destination is critically contingent upon the identifiability of 
(^, h) as a solution to Problem (Pi) by the criterion in Definition 1. 

As with our prior results, herein we shall consider /C = Pi x X >2 to be a separable cone in x with 
^ G Pi C and ^ G P 2 C In many applications, it is typical to have Pi (respectively P 2 ) to be a 
low-dimensional subset of R’^ (respectively R’^). An important interpretation of ^ G Pi is that g can be considered 
as a coded representation (of ambient dimension m) of the unknown whose actual dimension is much smaller (equal 
to that of the Hausdorff dimension of Pi). However, this coded representation is invoked to help strengthen the 
identifiability properties of the problem, rather than to help with error correction in the presence of observation 
noise. This is a somewhat different interpretation of coding than the classical goal of correcting noise induced errors. 
Operationally, both interpretations rely on redundancy albeit for different purposes. 


A. Repetition Coding 

Suppose that difew of the relays in Fig. 2 cooperatively decide to maintain the same relative transmission powers 
and phases and this decision is communicated to the destination for use as side information for channel estimation. 
A physical motivation for allowing few rather than all relays to cooperate, may be attributed to infeasibility of such 
coordination for widely separated relays and/or for a large number of relays. Mathematically, this translates to the 
use of a smaller feasible set JC by virtue of a stricter structural restriction on g via the set Pi that is smaller than R’^. 
Specifically, let A' C {/i, / 2 , • • •, h} ^ {1? 2,..., m} be the index set of propagation delays corresponding to the 
cooperating relays. Assuming that the only structure on ^ G Pi is that induced by the cooperating relays, we define 
the repetition coded domain as Pi = /Ci(A',m) where (A',m) act as parameters to the following parametrized 
definition of a family of cones for any dimension d>2 and any index set A C {1, 2,..., d}. 


ICi{A,d) ^ e 


^^(l) ^ 0, w{d) ^ 0, 3 c € M \ {0}, such that w{A) = cl|. 


(16) 


In other words, the set Pi imposes the structure that each of its member vectors has the same non-zero value 
c on the index subset A'; the value of c being allowed to vary across members. Even for a large number of 
cooperating relays (as measured by the cardinality |A'|), Corollaries 3 and 4 stated below imply unidentifiability 
results for h (unstructured channel) and h G /Co(A",n) (canonical-sparse channel), respectively. Corollary 3 
makes a statement for the repetition coded domain /Ci(A',m) that is analogous to the conclusions drawn by 
Corollary 1 for the canonical-sparse domain /Co(A,m). The main technical difference between Corollaries 3 and 1 
is that they require different constructions for a candidate unidentifiable input signal since the feasible domain 
JC is different across these results. We defer further technical details to the proof of Theorem 3 in Appendix C. 
Corollary 4 can be interpreted as extending the unidentifiability results of Theorem 2 from the canonical-sparse 
product domain /Co(Ai,m) x /Co(A 2 ,n) to the mixed product of repetition coded and canonical-sparse domains 
/Ci(A',m) X /Co(A",n). Roughly speaking, because we only consider separable domains JC, one can simply fuse 
the different constructions of a candidate adversarial signal x* G /Ci(A',m) from Corollary 3 and y* G /Co(A",n) 
from Corollary 1 to produce a candidate unidentifiable pair G /Ci(A',m) x /Co(A",n) for Corollary 4. This 

ability to fuse individual adversarial signal constructions on sets Pi and P 2 to produce candidate unidentifiable 
signal pairs within the Cartesian product domain /C = Pi x P 2 underlies the simplicity of arguments in the proof 
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of Theorem 3. At present, we cannot foresee any adaptations to the constructions used in the proof of Theorem 3 
that would allow us to make analogous statements for non-separable domains /C in a generalizable way. 

Corollary 3. Let m > 3 and n > 2 be arbitrary integers. For any given index set ^ ^ A' C {2, 3,..., m — 1}, 
let 1C = /Ci(A',m) x and define p = |A'IJ(A' — 1)|. Then there exists a set C /C/Mr such that every 
{d: h) ^ JC is unidentifiable. If A' ~ 1) ^ 0 of dimension at least {m-\- n — p), otherwise is of 

dimension at least {m + n — p + 1). 

Corollary 4. Let m > 3 and n > 5 be arbitrary integers. For any given index sets 0y^A' C — 1} and 

0 A" C {3,4,..., n — 2}, let 1C — /Ci(A', m) x /Co(A", n) and let pi = |A' ~ 1)1 P2 — 1^^^ ~ 1)1- 

Then there exists a set C /C/ ITr such that every (g^h) G is unidentifiable. If A' ~ 1) ^ 
dimension at least {m -\- n — pi — P 2 ), otherwise is of dimension at least (m + n + 1 — pi — ^ 2 )- 

We note that in contrast to the results for the canonical-sparse domain /Co(A,m) in Corollary 1, the dimension 
of unidentifiable subset for the repetition coded domain /Ci(A',m) as given by Corollary 3 depends on whether 
A' ~ 1) is non-empty, i.e. if atleast one pair of adjacent indices for g E /Ci(A', m) are repetition coded and 
hence belong to A'. Corollary 4 essentially says that even with a canonical-sparse prior structure on the channel 
vector h and nearly full cooperation between intermediate relays (|A'| = 0(m)) via repetition coding on g, there 
may be non-zero dimensional unidentifiable signal subsets in the feasible domain 1C. We recall that the index subsets 
A' and A" in Corollary 4 are known at the receiver and the unidentifiable signal set exists despite the availability of 
this side information. Thus, the repetition coding based system architecture laid out in this section is bound to fail 
in practice although it may be an intuitive choice. An important implication of Corollary 4 is that it is necessary 
to have A'fj{l,m}y^0, which is equivalent to the requirement that either the first relay or the last relay (when 
ordered by increasing propagation delays) must be present in the cooperating subset of relays. 


B. Partially Cooperative Coding 

The results in Section IV-A can be generalized to other codes. Consider the system described in Section IV-A, 
except that the cooperating relays decide on a different structural restriction for g ^ Vi C As before, let 
A' C {1,2,... ,m} denote the index set of propagation delays corresponding to the cooperating relays. We let 
= lCb{A'^m) denote a partially cooperative coded domain where the vector b E I denotes the cooperative 
code (alternatively, b could also be interpreted as a power-delay-profile across cooperating relays) and (A',m) are 
specific parameters to the following parametrized definition of a family of cones for any dimension d>2 and any 
index set A C {1,2,..., d}. 


ICbiA, rf) = e 


tu(l) ^ 0, w{d) ^ 0,3c € M \ {0}, such that w{A) — c6|. 


(17) 


As a visual example. Fig. 3 shows an arbitrary vector s in the partially cooperative coded domain lCb{A^d) = 
lCb{{3, 4, 7, 8, 9,12}, 14) for the code vector b = (0.5,0.835, —0.3, —0.5, —0.835, —0.15)^. It is important to note 
that the code vector b E I is known at the destination so that the decoder has explicit knowledge of the 
set Vi = lCb{A'^m). Furthermore, the partially cooperative coded domain lCb{A^d) unifies and generalizes both 
the repetition coded domain /Ci(A,(i) and the canonical-sparse domain lCo{A^d). To see this, we first note that 
6 = 0 implies lCb{A^ d) = /Co(A, d) by definition, and that the definitions (16) and (17) are equivalent for 6 = 1. 
Secondly, if the code vector 6 E ]RI^iUA 2 | admits a partition of the form 6 (Ai)^, 6 (A 2 ) ^ = [ 0 ^, 1 ^] for 


some disjoint index subsets Ai, A 2 C {1,2,... ,(i}, then the definitions in (17), (16) and (6) would imply that 
U ^2, G?) = /Co(Ai, (i) fj /Ci(A 2 , d). This can be checked directly by substitution. 

Corollary 5 and Theorem 3 stated below imply unidentifiability results for certain joint properties of the (A', 6) 
pair, i.e. for certain joint configurations of the cooperating relay index set and the cooperative code employed. The 
unidentifiability results may hold for both 6- E (unstructured channel) and h E lCt,{A"^n) (subspace-sparse 
channel), even for a large number of cooperating relays (as measured by the cardinality |A'|). In what follows, 
A* + 1 denotes the Minkowski sum of the sets {1} and A*. 


Corollary 5. Let m > 3 and n > 2 be arbitrary integers. For any given index set ^ A' C {2, 3,..., m — 1}, let 
6 0, /C = lCb{A' ^ m) X MF and define p = | A' ~ 1) I- — 1) = 0, then there exists a set G^ C /C/ Mr 

of dimension at least {m + n — p + 1), such that every (^, h) E G^ is unidentifiable. Otherwise, if A' ~ 1) 7^ ^ 
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index, I 


Fig. 3. An arbitrary vector s G /Cb(A,d) with A = {3,4,7,8,9,12}, d = 14 and b = (0.5,0.835,—0.3,—0.5,—0.835,—0.15)^. The 
specific vector s in the left plot satisfies s(A) = cb with c = — 1, verifying (17) with code vector b as in the right plot. Every s G /Cb(A, d) 
is collinear with b on the index set A (indicated by red stems) with A representing the identities of the cooperating relays. Heights of the 
black dashed stems, indicating values on A^^, represent contributions from non-cooperating relays and can vary across different vectors in 
/Cb(A,d). 


then denote 0 ^ A* C {1, 2,..., |A'|}, the index subset such that Mg G /C6(A', m), ~ collinear 

with If and b{A^ + 1) are collinear, then there exists a set C /C/Mr of dimension at least 

{m + n — p), such that every (g^h) G is unidentifiable. 

Corollary 5 makes a statement for the partially cooperative coded domain /C 5 (A',m) that is essentially analogous 
to the statement of Corollary 3 for the repetition coded domain /Ci(A',m), the only technical distinction being a 
more elaborate specification of the unidentifiable configurations of the pair (A', b) in the former result. Some of the 
technical conditions on the pair (A', 6), for generating unidentifiable configurations, are automatically satisfied when 
5=1 and hence we chose to present this simpler instantiation of Corollary 5 in the form of Corollary 3 under 
the rubric of repetition coding in the last subsection. In particular, both domains ICtiA'^m) and /Ci(A',m) share 
the same construction for a candidate unidentifiable input signal, provided the premise of Corollary 5 is satisfied. 
Furthermore, Corollaries 1, 3 and 5 motivate the following categorization of the (A, 5) pairs in preparation for the 
statement of Theorem 3. 

Definition 2. Let d > 3 be an arbitrary integer, 0y^AC {2,3,...,(i—1} be a given index set and 5 G 5^ a 
code vector. Let p = |AIJ(A — 1)|. The following mutually exclusive (but not exhaustive) categories for the pair 
(A, b) are defined. 

0) If d > 5, 0 A C {3,4,..., (i — 2} and 6 = 0 then (A, 6) is of type 0. 

1) Let 6 0 and A P|(A — 1) 0. is collinear with 1 G then (A, 6) is of type LIfb is not collinear with 

1 G then let 0 A^, C {1, 2,..., |A|} denote the index subset such that Mw G JCh{A^ d), iu(A P|(A — 1)) 

is collinear with 6(A*). Ifb{Af) and 6(A* + 1) are collinear, then (A, 6) is of type 1. 

2) Ifby^O and A P|(A — 1) = 0 then (A, 6) is of type 2. 

For a visual example, we note that Fig. 1 represents a (A, 6) pair of type 0 while Fig. 3 shows a type 1 pair with 
A,^ = {1,3,4} (A* as defined above in Definition 2). 

Theorem 3. Let m^n > 3 be arbitrary integers. For any given index sets 0 A' C {2, 3,..., m — 1} and 
0 7 ^ A" C {2,3,..., n — 1}, let 1C — lCh{^'•> '^) x lCb'{A'\ n) and let p = |A' “ 1)1 P' — W' “ 1)1- 

If the pairs (A', 6) and (A", 6') are respectively of types t and f for G {0,1,2}, then there exists a set 
G^ C /C/ Mr of dimension at least {m + n — 1 — p — p' + t + f) such that every (^, h) G G^ is unidentifiable. 

Proof: Appendix C. ■ 

Theorem 3 stated above utilizes the construction of candidate adversarial inputs over different domains, as developed 
in Corollaries 1 and 5, and fuses them together to generate candidate unidentifiable signal pairs by exploiting the 
separability of the feasible domain 1C. It is straightforward to check that Theorem 2 and Corollaries 1, 3, 4, and 5 
can all be derived from Theorem 3 by suitable choices of the parameter pairs (A', 6) and (A", 6'), e.g. setting 6 = 1 
and 6' = 0 gives back Corollary 4. In essence, barring Corollary 2 which requires special assumptions. Theorem 3 
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encompasses all of the unidentifiability results for constrained blind deconvolution developed in this part of the 
paper. 

The main message of Theorem 3 is that the coded subspaces represented by the vectors b and b' are critical to 
the identifiability of the blind linear deconvolution problem (Pi). If these coded subspaces are solely determined by 
the application or natural system configuration (like the multi-hop channel estimation example above) then sparsity 
may not be sufficient to guarantee identifiability for the blind deconvolution problem. A certain amount of design 
freedom on both of the coded subspaces (represented by vectors b and b') is necessary to guarantee identifiability 
under blind deconvolution (like in [5], [14], [20]). 

As yet another important special case for wireless communications, we interpret the vector 6 as a power- 
delay profile and assume it to be geometrically decaying over some contiguous index subset A, i.e. b satisfies 
b{j) = rb{j — 1) for some 0 < |r| < 1 and all j G {2,3,..., |A|}. Then we get the geometrically decaying cone 
lCiy{K^d) of all (i-dimensional real vectors that are geometrically decaying on the index subset A and it admits the 
following unidentifiability result (in the statement below, /C 5 (A 2 ,n) represents a geometrically decaying cone). 

Corollary 6. Let m > 3 and n > 5 be arbitrary integers. For any given index subset 0 ^ Ai C {3,4, ..., m — 2} 
and any contiguous index subset 0 ^ A 2 C {2,3,..., n — 1}, let fC — /Co(Ai,m) x /C 5 (A 2 ,n) and define 
Pi = |Ai U(Ai — 1)|, p 2 = IA 2 I + 1. Then there exists a set Q C /C/Mr such that every {g^h) G Q is 
unidentifiable. ^ IA 2 I = 1 then Q is of dimension at least {m + n+1 — pi — P 2 ), otherwise Q is of dimension at 
least {m-\- n — Pi — P 2 ). 

We omit the proof of Corollary 6 here as it follows from specializing the proof of Theorem 3. A complete proof 
of Corollary 6 appears in our earlier work [2]. The occurrence of a geometrically decaying power-delay profile is 
common in underwater acoustic communication channels [44] and more generally in wide-band communication 
channels [45], [46]. 


V. Conclusions 

Blind deconvolution is an important non-linear inverse problem commonly encountered in signal processing 
applications. Natively, blind deconvolution is ill-posed from the viewpoint of signal identifiability and it is assumed 
that application specific additional constraints (like sparsity) would suffice to guarantee identifiability for this inverse 
problem. In the current work, we showed that (somewhat surprisingly) sparsity in the canonical basis is insufficient 
to guarantee identifiability under fairly generic assumptions on the support of the sparse signal. Specifically, we 
explicitly demonstrate a form of rotational ambiguity that holds true for sparse vectors as well. Our approach builds 
on the lifting technique from optimization to reformulate blind deconvolution into a rank one matrix recovery 
problem, and analyzes the rank two null space of the resultant linear operator. While this approach is philosophically 
applicable to other bilinear inverse problems (like dictionary learning), it is the simplicity of the convolution operator 
in the lifted domain that makes our analysis tractable. We developed scaling laws quantifying the ill-posedness of 
canonical-sparse blind deconvolution by bounding the dimension of the unidentifiable sparse signal set. When applied 
to a second-hop sparse channel estimation problem, our methods also revealed the insufficiency of side information 
involving repetition coding or geometrically decaying signals towards guaranteeing identifiability under blind linear 
deconvolution. To establish our results, we developed a measure theoretically tight, partially parametric and partially 
recursive characterization of the rank two null space of the linear convolution map. This result is a precursor to 
non-randomized code design strategies for guaranteeing signal identifiability under the bilinear observation model of 
linear convolution. The design of such codes is a topic of ongoing research. 
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Appendix A 
Proof of Theorem 2 


The proof relies on constructing a generative model for vectors x G /Co(Ai,m) and y G /Co(A 2 ,n) such 
that (x^y) is unidentifiable within JC = /Co(Ai,m) x /Co(A 2 ,n). Let u G /Co(Ai |J(Ai — 1),m — 1), v G 
/Co(A 2 U (^2 — l),n — 1), and (0,0) C Ga ^ [0,27r)^ be chosen arbitrarily, where 


Qa = |(/?,7) e [0,27r)^ 


^,7 ^ {/7r/2 \l eZ}, P 


7^p7r \leZ}y 


(18) 


We use the 4-tuple (rx, u, 0,0) to generate the vector pair (x, y) using 


X = 

u 0 

0 —u 

cos 9 
sin 9 

, y = 

0 V 

V 0 

sin0 

— COS0 

, (19) 

and the vector pair (x', y') using 

/ 

u 0 

COS0 

/ 

0 V 

sin0 

(20) 

X = 

0 —u 

sin0 

, y = 

V 0 

— COS0 


By assumption, 0 0 {i6(l),rx(m — 1),u(l),u(n — 1)} and 0 0 {sin0,cos0,sin0,cos0}, so (19) implies 0 0 
{x(l), x(m), y(l), y(n)} and (20) implies 0 0 {x'(l), x'(m), y'(l), ^(n)}. We also have rx(Ai) = rx(Ai — 1) = 0 
and 'u(A 2 ) = v{K 2 — 1) = 0 by assumption, and 


x(Ai) = rx(Ai) COS0 — u{Ai — 1) sin0 = 0, y{A 2 ) = —u(A 2 ) cos0 + u(A 2 — 1) sin0 = 0, (21a) 

x'(Ai) = u{Ai) COS0 — u{Ai — 1) sin0 = 0, y{A 2 ) — —v{A 2 ) cos0 + v{A 2 — 1) sin0 = 0, (21b) 


where (21a) follows from (19) and (21b) follows from (20). Therefore, x,x' G /Co(Ai,m) and y^y' G /Co(A 2 ,n), 
implying that (x,y), {x'^y') G 1C. Setting X = xy^ — x\y')^, we have 


X = 


+ 


u 

0 ■ 


COS 9 

r 

0 

—u 


sin 9 

[ 

u 

0 ■ 


COS0 

r 

_0 

—u 


sin0 

[ 


[sin 0 — cos 0] 

[— sin 0 cos 0] 


0 

T 

0 

.,T 


= sin(0 — 0) 


u 0 


0 

o 

1 

1 _ 


O 


V 


( 22 ) 


and hence, xy^ — x'(y')^ = X G A/^(^,2) from Lemma 1. Therefore, the pairs (x,y) and (x',y') produce the 
same convolved output and are indistinguishable under the linear convolution map. Since x and x' are linearly 
independent by the choice (j) — 9 ^ {Itt \ I G Z}, (x, y) C JC is unidentifiable by Definition 1 with (x', y') as the 
certificate of unidentifiability. 

All that remains is to lower bound the dimension of a set ^ C /C of pairs (x,y) that can be shown to be 
unidentifiable for blind linear deconvolution using the above mentioned construction. Let Qy denote the Cartesian 
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product set /Co(Ai lJ(Ai — l),m — 1) x /Co(A 2 U (^2 ~ 1 )?'^ ~ !)• We construct Q as follows. A vector pair 
{x^y) ^ Q if and only if all of the following conditions are satisfied. 

(Al) (x,y) G X is generated from a 4-tuple {u^v,9,(p) G x x [0,27r)^ using (19). 

(A2) (0,0) G^a. 

(A3) (u, r;) G Gy. 

From the arguments in the previous paragraph, we have 0 ^ C /Co(Ai,m) x /Co(A 2 ,n) = 1C and furthermore, 

every G ^ is unidentifiable within 1C. Let us consider the sets 

G'l = {('*^,0) ^ X [0,27r) I (A2) and (A3) are true}, 

^2(0) = {('^7 0) ^ X [0,27r) I given 0 G [0, 27r) \ {/7r/2 | I G Z}, (A2) and (A3) are true}. ^ ^ 

By the definition in (6), /Co(Ai U(Ai — 1), m — 1) is a (m — 1 — |Ai U(Ai — 1)|) = (m — 1 — pi) dimensional 
Borel subset of and /Co(A 2 U (^2 — 1), n — 1) is a (n — 1 — IA 2 U (^2 ~ 1 ) 1 ) = (n — 1 — P 2 ) dimensional 

Borel subset of Further, [0,27r) \ {/7r/2 | / G Z} is a one dimensional Borel subset of R and given a value 

of 0 G [0, 27r), [0, 27r) \ {/7r/2, 0 + /tt | / G Z} is also a one dimensional Borel subset of R. Hence, (23), (A2) and 
(A3) imply that 

1) is a (m — 1 — pi) + 1 = (m — pi) dimensional Borel subset of R’^, 

2) given 0 G [0, 27r) \ {/7r/2 | I G Z}, G 2 {^) is a (n — 1 — P 2 ) + 1 — {n — P 2 ) dimensional Borel subset of R’^, 

3) Ga is a two dimensional Borel subset of R^, and 

4) Gy X ^A is a (m — 1 — pi) + (n — 1 — P 2 ) + 2 = (m + n — pi — P 2 ) dimensional Borel subset of 
Next, we compute the dimension of the set G> Consider a factorization of G into 

Gi = {x^ G R"^ I {x^,y^) G G] and G 2 {x^) = {y^ G R"" | (x*,?/*) G G], (24) 

i.e. G = {(x*,y*) I Xx. G 01, C 02(x*)}. Equations (24) and (19) imply that every x G 0i is generated as the 
result of a uniformly continuous map from (it, 0) G Gi (with a non-singular Jacobian matrix). From Lemma 2, the 
quotient set Qr^{x^m) is finite for every x G 0i and therefore each x G 0i can be generated by at most a finite 
number of elements (it, 0) G G'l using (19). Given some x G 0i, let 0 C [0, 27r) be a set such that for every 0 G 0, 
there exists (it,0) G Gi generating x using (19). Clearly, |0| < (2m — 2) from Lemma 2. Then, (24) and (19) 
imply that every y G 02(x) is generated as the result of a uniformly continuous map from ( 1 ?, 0) G U6>g© ^ 2 (^) 
(with a non-singular Jacobian matrix). Using Lemma 2, the quotient set Q^(y, n) is finite for every y G 02 (x) and 
therefore each y G 02(x) can be generated by at most a finite number of elements ( 1 ?, 0) G U6>g© ^ 2 (^) using (19). 
Since 0 is a finite set, the above arguments and (24) imply that every element (x, y) G 0 is generated by at most a 
finite number of 4-tuples (i6,i?,0,0) G 0v x 0 a, using the uniformly continuous maps in (19) (with non-singular 
Jacobian matrices). Since 0v x 0 a is a (m + n — pi — P 2 ) dimensional Borel subset of R^+^, we get 0 as a Borel 
subset of R^+^ of dimension {m + n — pi — P 2 )> 

We let 0* = 0/ IdR, i.e. G^ is the quotient set of G w.r.t. the equivalence relation ITr. Clearly, G^ ^ 1C j ITr. Since 
each element (x, y) G 0/ ITr is a representative for a one dimensional Borel subset {(ax, ^y )} C 0 C R^ x R’^, 
the dimension of 0/ ITr is one less than the dimension of 0. Hence, 0* is a (m + n — 1 — pi — P 2 ) dimensional 
set. 


Appendix B 
Proof of Corollary 2 

This proof follows the template of the proof for Theorem 2 in Appendix A with adjustments for the change 
of feasible set 1C and borrows from the proof of Theorem 1 in Part I of the paper to show almost everywhere 
unidentifiability. 

Let u G /Co(AU(A — l),m — 1) and y G R’^ be arbitrarily chosen. Since n > 4 is an even integer, invoking 
Lemma 2 for y yields a vector v G R^“^ and a scalar 0 G [0,27r) such that the second relationship in (19) is 
satisfied. Next, we select 0 G [0, 27r) \ {/7r/2 | / G Z} and use the parameters (u, 0) in (19) to generate the vector 
X G R’^. We reuse the parameter 4-tuple ( 16 , 1 ;, 0,0) in (20) to generate the pair (x',y') G R’^ x R’^. Note that 
our construction for the 4-tuple (u, r;, 0, 0) here is different from that in Appendix A, even if we set A 2 = 0 in 
Theorem 2. Nonetheless, using (22) and reasoning on the same lines as in Appendix A leads to the conclusion that 
(x,y) G /C is unidentifiable with (x',y') as the certificate of unidentifiability. 
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For notational brevity, let us define 

7^^/Co(a|J(A-1) ^ 777/ ^^ ^ 277^ ^ "1^/77^2 | / ^ ^^ . (^25^ 

To show that given any {u^9) G %, almost every choice of y E yields an unidentifiable pair {x^y) E /C, it 
suffices to show that (j) 0 {17t/2^6 + /77 | / E Z} holds almost everywhere w.r.t. the measure over y. From (19), 
the dependence of y is uniformly continuous on the n real numbers (0,17(1), 17(2),... ^v{n — 1)) and these n 
real numbers completely parametrize y. Further, since the measure over y is absolutely continuous w.r.t. the n 
dimensional Lebesgue measure, it is possible to choose a measure over (j) that is absolutely continuous w.r.t. the 
one dimensional Lebesgue measure. Since {/77/2, 0 + /77 | / E Z} is a zero dimensional set (given the value of 9), 
(j) 0 {/77/2,0 + /77 I I E Z} is true almost everywhere w.r.t. the measure over (j) and therefore also w.r.t. the measure 
over y. 

Let T-L' C /Co(A, m) denote the image of T-L under the (rx, 9) ^ x map in (19). To complete the proof, we need to 
bound the dimension of the set H'. Analogous to Appendix A, we have /Co(A U(A — 1), ttT/ — 1) as a {m — 1 — p) 
dimensional Borel subset of and therefore, is a (ttT/ — 1 — p) + 1 = (ttT/ — p) dimensional Borel subset of 

W^. For any x E H', Lemma 2 implies the finiteness of the quotient set Qr^{x^m) and therefore each x ^ T-L' 
can be generated by at most a finite number of elements (xx,0) E T-L using (19). Since (19) represents uniformly 
continuous maps with non-singular Jacobian matrices, we have T-L' as a (ttT/ — p) dimensional Borel subset of W^. 
The proof is complete since for any x E T-L', {x^y) E /C is unidentifiable almost everywhere w.r.t. any measure over 
y that is absolutely continuous w.r.t. the n dimensional Lebesgue measure. 


Appendix C 
Proof of Theorem 3 

For notational convenience, we denote g hy x, h by y, by and /C5/(A",77/) by V 2 {n). The 

high level architecture of this proof is similar to the proof of Theorem 2 in Appendix A. However, the generative 
model for unidentifiable pairs (x,p) E /C is substantially different, owing to the change in the feasible set /C. 

We shall use (19) to construct unidentifiable pairs (x,p) in the separable feasible set JC = Vi{m) x V 2 {n). 
Since the generative parameters for x and y in (19) are disjoint, an analogous construction also works for the 
feasible set K' — V2{m) x Vi{n) to produce unidentifiability results. For example. Corollary 4 asserts that for 
/C = /Cl(A', 777/) X /Co(A", 77/) and A'P|(A' — 1) 0, there exists an unidentifiable subset of /C/Mr of dimension 

(ttT/ + 77/ — Pi — P 2 ). Our reasoning based on disjoint generation would imply that for 1C' = /Co(A', m) x /Ci(A", n) 
and A" n(A^^ — 1) 0, there exists an unidentifiable subset of fC'/ ITr of dimension {m + n — pi — P 2 ). Thus, it 

suffices to consider unordered pairs (t, t') E {0,1, 2}^ for the purpose of this proof. Note that Theorem 2 is exactly 
equivalent to the case of (t, t') = (0, 0). In the rest of the proof, we complete the treatment for the remaining cases 
(t,Oe{0,l,2}E{(0,0)}. 

Let us choose the 4-tuple of parameters (m, u, 6, cp) G R™ ^ x ^ x [0,27r)^ such that we generate (x, y) G 
1C — lCi){A'^m) X /C5/(A",77/) from (19) and {x'^y') E 1C from (20). For the construction to be consistent, u must 
satisfy 


x(A') = '^(A^) cos 9 — u(^A' — l) sin 9 = ci(0,0)5, (26a) 

x'(^A') = 'r^(A') COS0 — u(^A' — l) sin0 = C 2 (0,0)6, (26b) 


for some scalars ci(0,0),C2(0,0) E M \ {0} that could depend on {9^(1)), and 0 0 {u{l)^u{m — 1)}. Solving (26) 
for u{A') and u{A' — 1), we get 



Cl (0, 0 ) sin 0 — C 2 (0, 0 ) sin 0' 
sin(</) - 6) 

ci{6,4>) cos 4> — C 2 { 0 ,(p) cos 9 
sin(^ — 9) 


b, 

\b. 


Analogously, consistency also requires v to satisfy 


(27a) 

(27b) 


y (A") = v[A" — l) sin^ — ^(A") coscp = c'i{9, (p)b', (28a) 

y'[A") — v(^A" — l) sin0 — 'u(A'') cos9 — C2{9, (t>)b\ (28b) 
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for some scalars c'^( 6 ^, 0), C 2 ( 6 ^, 0) G M \ {0} that could depend on (0, 0), and 0 0 {v{l)^v{n 
for v{h") and v{h" — 1 ), we get 


V 


t,(A») = ( 
(A" - 1) = ( 


C|( 0 , (j)) sin 0 — C 2 ( 0 , (j)) sin(/)\ , 
sin((/) — 0 ) y ’ 

C|(0, 0) COS0 — C2(0, 0) COS(/)\ , 

sin((/) — 9) J 


1)}. Solving (28) 


(29a) 

(29b) 


Given a value of (0, (p) G (Ga is defined in (18)), if there exist solutions to u satisfying (27) and v satisfying (29), 
then the construction of (x,y), {x'^y') G /C is consistent. Invoking (22) and using the same line of reasoning as in 
Appendix A implies that {x^y) G /C is unidentifiable with {x^y') as the certificate of unidentifiability. Furthermore, 
the consistency condition (27) for u is independent of the consistency condition (29) for v for a given (0,0) G Qa. 
Thus, it suffices to analyze the cases t G {0,1,2} for different values of (0, 0) and the same analysis would apply 
to the cases t' G {0,1,2} owing to the similarity between (27) and (29). 

If (A', b) is of type t = 0, the analysis in Appendix A can be reused. In particular, we have 6 = 0 and (27) 
reduces to the requirement u(A' 1J(^^ ~ 1)) — Clearly, any u G /Co(A' [J{A' — 1), m — 1) is a valid assignment 
and /Co(A' |J(^^ — 1), m — 1) is a (m — 1 — p + t) dimensional Borel subset of for t = 0. 

If (A', 6) is of type t = 2, then A' — 1) = 0 and there are no common variables between the vectors u{A') and 
u{A' — 1). Thus, any 0 0 {0 + stt | 5 G Z} will yield a consistent assignment for u(^A' U(A' — 1)) through (27) with 
arbitrary constants ci(0,0) and C 2 ( 0 ,0) satisfying ci(0,0) sin0 7 ^ C 2 ( 0 , 0) sin0 and ci(0,0) cos0 7 ^ C 2 ( 0 ,0) cos0. 
Clearly, for (0,0) G Ga, there is a two dimensional set of valid choices for the pair (ci(0,0),C2(0,0)) and 
therefore u{A') and rx(A' — 1) lie (independently of each other) on a one dimensional subspace collinear with 
6 G implying that 16 (A'IJ(A' — 1)) lies on a two dimensional subspace of W. Since u{j) is unconstrained 
Vj 0 A'|J(A' — 1), the set of consistent assignments of u is a Borel subset of of dimension at least 

2 + (m — 1— p) = {m — 1 — p + t). 

Finally, let us analyze the case when (A', 6 ) is of type t = 1. We have A' f)(A' — 1 ) 7 ^ 0, so u(A' P(A' — 1 )) 
must receive consistent assignments from (27a) and (27b). If 6 is collinear with 1, then (27) is consistent if and 
only if sin (0 — 0 ) 7 ^ 0 and 


ci( 0 , 0 ) sin 0 — C 2 ( 0 , 0 ) sin 0 = ci( 0 , 0 ) cos 0 — C 2 ( 0 , 0 ) cos 0 7 ^ 0 . 


(30) 


If 6 is not collinear with 1 , then there exists an index subset 0 7 ^ A^ C {1,2,..., |A'|} as characterized in 
Definition 2. Further, (27a) and (27b) respectively imply (31a) and (31b) below. 


«(A'n(A'-i)) 

«(A'n(A'-i)) 


/ci{e,(()) sin(f) - 02(9, (p) sine\ 

I- 

f ci{ 9 ,(j)) COS(p - 02(9, (p) COS 9 \ 

[ --ji>(A. + l) 


(31a) 

(31b) 


Since (A', 6 ) is of type t = 1, b{A^) and 6 (A* + 1) are collinear and there exists a constant r G M \ {0} such that 
6 (A* + 1) = r 6 (A*). The consistency requirement in (31) then reduces to requiring sin(0 — 0) 7 ^ 0 and 


ci( 0 , 0 ) sin 0 — C 2 ( 0 , 0 ) sin 0 = r • ci( 0 , 0 ) cos 0 — r • C 2 ( 0 , 0 ) cos 0 7 ^ 0 . 


(32) 


Note that (32) subsumes (30) (on setting r = 1), and further implies 

C 2 ( 0 , 0 ) sin0 —rcos0 sin(0 — S) 
Cl (0,0) sin0 —rcos0 sin(0 — S) 

for B — arctan(r). We define the set — [O’ 2^)^ 


= {(^> 7 ) e [0,27r)^ 


0,7 0 {17t/2^ / 7 r + S|/GZ},0 — 70 {In | / G Z}|. 


(33) 


(34) 


and notice that for any (0,0) G 0a(^)’ sin(0 — 0) 0 0 and (33) implies that the set of valid choices for the pair 
(ci(0,0), C 2 ( 0 , 0)) is one dimensional. Therefore, it follows from (27) and (31) that rx(A'lJ(A' — 1 )) lies on a one 
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dimensional subspace of W and is given by 


ru(A'\(A'-l))] 


■ ' 

«(A'n(A'-i)) 

[u((A'-l)\A')J 

= c(6',0) 

6 (A*) 


for some scalar c(0,0) G IR\ {0}. In the special case when b is collinear with 1, (35) simplifies to u(A' 1J(^^ ~ 1)) 
being collinear with 1. Since u{j) is unconstrained Vj 0 A' IJ(A' — 1), the set of consistent assignments of u is a 
Borel subset of of dimension at least 1 + (m — 1 — p) = (m — 1 — p + f), for any (0, 0) G 

For notational brevity, we define the following. Let C (respectively Q 2 ^ denote the set of all 

vectors u (respectively v) that satisfy the consistency condition (27) and 0 0 {u{l)^u{m — 1)} (respectively (29) 
and 0 0 {i;(l),i;(n — 1)}) for a given value of (0,0). Moreover, let Qs C M? denote the largest set of allowed 
values for the pair (0,0) depending on the type pair (f, t') of the pairs (A', b) and (A", b'). In particular, 

1) G { 0 , 2 } implies Gs ^ Ga. 

2 ) f = 1 implies G ^ ^a(^) ^ and 

3) t' — 1 implies G ^ ^ where r' depends on (A", b') in the same way as r depends on (A', b). 

Finally, let G' ^ 1C denote the image of the set Gi x G 2 x Gs under the map (u^v^O^c/)) i-G (x,p) in (19). 
Clearly, ^ 0 and every (x*,p*) G G' is unidentifiable within 1C if (27) and (29) can be simultaneously 
satisfied with 0 0 {u{l)^u{m — 1),'u(l),'u(n — I),sin0,cos0,sin0,cos0}. Our analysis indeed guarantees this 
using Definition 2 and careful descriptions for the sets Gi, G 2 , Gs and the conic family in (17). 

To finish the proof, we need to compute the dimension of G'- Our analysis for the cases t G {0,1,2} establishes 
that is a (m — 1 — p + f) dimensional Borel set and (by analogy) ^2 is a (n — 1 — p' + 1 ') dimensional Borel 
set. Clearly, Gx{r) is a two dimensional Borel set for any r G M. Therefore, x ^2 x ^3 is a (m — 1 — p + f) + 
(n — 1 — p' + f') + 2 = {m -\- n — p — p' -\- 1 -\- 1 ') dimensional Borel subset of Since G' is the concatenated 

image of the uniformly continuous maps in (19) (with non-singular Jacobian matrices), invoking Lemma 2 as in 
Appendix A implies that is a {m + n — p — p' + t + f) dimensional Borel set. Setting G^ C G' j Idn C /C/ Mr 
and invoking the argument in the last paragraph of Appendix A implies that the dimension of G^ is one less than 
the dimension of G' and completes the proof. 






